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ABSTRACT 

One key piece of information missing from high redshift galaxy surveys is the galaxies' 
cold gas contents. We present a new method to indirectly determine cold gas surface 
densities and integrated gas masses from galaxy star formation rates and to separate 
the atomic and molecular gas components. Our predicted molecular and total gas 
surface densities and integrated masses are in very good agreement with direct mea- 
surements quoted in the literature for low and high-z galaxies. We apply this method 
to predict the gas content for a sample of ~ 57000 galaxies in the COSMOS field at 
0.5 < z < 2.0, selected to have Iab < 24 mag. This approach allows us to investigate 
in detail the redshift evolution of galaxy cold and molecular gas content versus stellar 
mass and to provide fitting formulae for galaxy gas fractions. We find a clear trend 
between galaxy gas fraction, molecular gas fraction and stellar mass with redshift, 
suggesting that massive galaxies consume and/or expel their gas at higher redshift 
than less massive objects and have lower fractions of their gas in molecular form. The 
characteristic stellar mass separating gas- from stellar-dominated galaxies decreases 
with time. This indicates that massive galaxies reach a gas-poor state earlier than 
less massive objects. These trends can be considered to be another manifestation of 
downsizing in star formation activity. 
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1 INTRODUCTION 

In a pioneering work, iMadau et all (|l996t ) demonstrated 
that the star formation (SF) activity of the Universe peaks 
at a redshift z ~ 1 — 3 and decreases at lower redshifts to 
its present day value. In the following years, many studies 
confirmed the existence of a peak in the star formation rate 
(SFR) densi ty around these redshifts, and a clear dec line at 

2 < 1 (e.g.. lHopkinsll2004 [Hopkins fc Beacomll2006l ). This 
makes the epoch between z ~ 2 and z ~ 0.5, in which mas- 
sive galaxies stop forming the bulk of their stars and become 
passive, a crucial period for galaxy evolution. 

One of the features which has become clear in 
this period of galaxy evolution is that the most mas- 
sive objects (giant elliptical galaxies hosted in galaxy 
groups and clusters) formed their stars early, while 
less massive objects continued to form stars until the 
present (this trend is commonly called 'galaxy down- 
sizing'). This behavior has many different observa- 
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but its physical origin remains unclear 



|Fontanot et al.ll2009l i 

The star formation rate in galaxies is closely linked 
to the galaxy gas content. Observations have shown that 
SF in the Milky Way takes place i n dense, massive and 
cold giant molecul a r clouds (GMC; Solomon et al.l 1 19871 : 
iBolatto et all 120081 : iMcKee fc Ostrikerl 120071 1. This makes 
the SF ti ghtly con n ected to the molecular and atomic gas 
available. ISchmidtl l| 19591 ) found a power-law relation be- 
tween the surface density of SFR and gas surface den- 
sity. This work has become the cornerstone of a wealth 
of studies relating SFR surface densities and cold gas and 
molecular hydrogen surface densities (e .g., lKennicutt| [l998b: 
iBigiel et aiTl2008l : ISchruba etHHl201lh . Both the observa- 
tions of GMCs, and the power-law relations between SFR 
and gas densities demonstrate that information about the 
gas and its partition into atomic and molecular hydrogen is 
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essential for a proper understanding of galaxy evolution and 
the build up of stellar mass. 

An essential landmark and significant step forward in 
the study of galaxy evolution at higher redshifts has been the 
development of large samples of distant galaxies with exten- 
sive multi-wavelength information (e.g., GOODS and COS- 
MOS ). These samples comprise imaging data from X-rays 
to the infrared, which allow us to derive numerous physical 
properties of galaxies, such as stellar masses and star for- 
mation rates. However, one of the key pieces of information 
missing in these studies is the cold gas content. Deriving 
the amount of cold gas present as a function of stellar mass 
for large representative samples would allow us to better 
understand the history of gas consumption and the overall 
phenomenon of galaxy downsizing. 

Information about the gas content of galaxies also 
provides important constraints on theoretical models of 
galaxy formation. Gas fractions help to break the degen- 
eracies in different physical mechanisms that are included 
in these models, namely star formation and stellar feed- 
back (Caviglia & Somerville, in prep). Moreover, cosmolog- 
ical galaxy formation models have begun to include recipes 
that allow a detailed tracking of atomic and molecular 
gas with redshift in the context of general galaxy prop- 
erties (e.g., gas m ass functions, gas properties as a func- 
tion o f stellar mass; Obrcschkow & Rawlings 200 3; iFu et al.l 
20ld; iGnedin fc Kravtsovl l201ll; iKrumholz fc Dekell l201lF 



Kuhlen et a l. 2011; Lagos et al. 2 011bl lal. Popping et al., in 



prep). Observational constraints on atomic and molecular 
gas content as a function of galaxy properties and cosmic 
time are crucial for developing and calibrating these new 
models. 

Direct observational measurements of the gas content 
in distant galaxies are currently available for only a very 
limited number of galaxies, and are likely biased towards 
the most gas-rich objects. Gas is usually detected using 
CO molecules as a tracer, so no direct information about 
the atomic gas content of high redshift galaxies is currently 
available. In addition, the conversion from CO to molec- 
ular hydrogen mass is notor iously uncertain an d may de- 
pend on galaxy properties (|Genzel et al.l |2010| . and refer- 
ences therein). In the coming decade, it is hoped that new 
facilities like ALMA (Atacama Large Millimeter array) and 
the SKA (Square Kilometer Array) will reveal the gas con- 
tent in representative samples of high-redshift galaxies. 

In the meantime, we can obtain indirect constraints 
on gas content by using the observational estimates of 
galaxy size and SFR and empirical correlations between SFR 
density and gas den sity. This approach ha s been used by 
lErb et al.l (^POSh and lMannucci et al.l |2009r i, who obtained 
estimates of gas masses for high redshift galaxies by invert- 
ing the Schmidt-Kennicutt relation (KS law) described in 
Kennicutt (1998b), which relates SFR surface density to the 
combined atomic and molecular hydrogen surface density. 
However, the indirect estimates of gas fraction that they 
obtain do not agree well with the direct measurements, as 
we later show. This could be because the direct estimates 
are biased high, or it could be that the "total gas" form 
of the KS law does not apply. It has been shown that the 
KS law breaks do wn at the lowest gas surface densities 
ijBigiel et al.l bOoH - ). This breakdown is ascribed to the in- 
ability of the cold gas at low surface densities to collapse 



gravitationally and form molecular clouds from which stars 
may originate. SFR surface density correlates in an almost 
linear fashion with the molecular gas surfa ce density without 
breaking down at lower s urface densities |Bigiel et al.ll2008l . 
l201ll;ISchru"ba" et al.ll201ll ). Using an inverted molecular-gas- 
based SF law therefore seems more appropriate to probe the 
cold gas content of high-redshift galaxies. Furthermore, in- 
verting a molecular-gas-based SF law allows for compari- 
son between directly and indirectly measured molecular gas 
masses of high-redshift galaxies. 

In this paper we present an improved method to derive 
the gas properties of distant galaxies from their star forma- 
tion rate densities. An integral part of our new method is a 
prescription to calculate the molecular fraction of cold gas 
based on the results by iBlitz fc Rosolowskvl (|200fj| ). These 
authors found that the ratio between the molecular and 
atomic hydrogen surface density (Rn 2 = £h 2 /£hi) in discs 
can be described empirically as a function of the hydrostatic 
mid- plane p ressure P m , driven by the stellar and gas den - 
sity (see also lWong fc Blitzll2002l ; IBlitz fc Rosolowskvll2004l l. 
This allows us to invert the much tighter and more physi- 
cal molecular gas based star fo r matio n relation derived for 
nearby galaxies by iBigiel et al.l (|2008T ) , instead of the tradi- 
tional KS relation. 

We apply our method to ~ 57000 galaxies from the 
COSMOS survey in the redshift range of 0.5 ^ z < 2.0 to 
study the structure and properties of gas and the galaxy 
evolution over this important cosmic epoch. COSMOS was 
chosen for our study because it has well established disc sizes 
based on high resolution (Hubble Space Telescope) imaging, 
crucial for our method, imaging running from the X-ray to 
radio wavelengths to derive stellar masses and SFR, and 
spans a wide range of redshift. 

The outline of this paper is as follows. In Section [2] 
we present our method to indirectly measure the molecular 
and the total cold gas content in galaxies. In Section [3] we 
apply this method to both low- and high-redshift data with 
known gas properties and calibrate our results against direct 
measures of cold and molecular gas. In Section [4] we de- 
scribe the application of the developed method to ~ 57, 000 
galaxies from the COSMOS survey in the redshift range of 
0.5 ^ 2 ^ 2.0, and the resulting gas properties with respect 
to other galaxy parameters. We discuss our findings in Sec- 
tion [5] and summarize our results in Section [S] Throughout 
the paper we apply a ACDM cosmology with Ho — 70 km 
s _1 , f2 matt cr = 0.28, and Qa = 0.72. We assu me a uni- 
versa l Chabrier stellar initial mass function (IMF: Ichabrierl 
120031 ) and where necessary convert all observations used to a 
Chabrier IMF. All presented cold and molecular gas surface 
densities and integrated masses include a correction of 1.36 
to account for Helium. 



2 METHOD 

2.1 Calculating gas surface densities 

We infer gas surface densities using a combinati on of an em- 
pirica l molecular star formation law (based on lBigiel et al.l 
I2008T ) and a prescription to calculate t he fraction of molec- 
ular hydrogen in cold gas (based on IBlitz fc Rosolowskvl 
l2006h . 
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Table 1. Summary of the parameters used in our indirect gas measure method. These parameters lead to the overall best fit results for 
integrated masses in Section [3] 



Parameter 


Description 


Value 


Reference 


A S F 


Normalization of the star formation law (Mq yr — 1 kpc -2 ) 


1.09 X 10~ 2 


Sec. 1331 


jV SF 


Index which sets star formation efficiency 


0.5 


Sec. 1331 


^crit 


Critical surface density which sets star formation efficiency (Mq pc — 2 ) 


100 


Sec. [531 


^0 


Normalization of the Blitz & Rosolowsky (2006) power-law 


2.35 x 10" 13 


Lerov et al. (2008) 


a 


Index of the Blitz & Rosolowsky (2006) power-law 


0.8 


Lerov et al. (20081 


la 


Mean ratio between gas and stellar vertical velocity dispersion in the stellar disc 


0.4 


Elmeereen (19931 


Xgas 


Scale radius of the gas disc, relative to the stellar disc 


1.5 


Sec. [531 



We have slightly adap ted the star formation law de- 
duced bv lBigiel et al.l ()2008l ) to allow for higher star forma- 
tion efficiencies in hig h gas surface densi ty re gions. This is 
based on the results of lDaddi et al.l |2010l ) and lGenzel et all 
who found the star-formation at high surface densi- 
ties tc^Jol_low_the KS - law (power- law slope of 1.4 versus 1.0 
for iBigiel et al.1 i|2008l )). The resulting equation is given by 

ESFH = IOaJpc-^ i 1 + fe) /h = Egas (1) 

where Esfr and E gas are the star formation and cold gas sur- 
face densities in Mq yr _1 kpc -2 and Mq pc -2 , respectively, 
^4sf is the normalization of the power law in Mq yr _1 kpc -2 , 
E cr it a criti cal surface density above which the star forma- 
tion follows iKennicutH (|l998bl ). Ng-p is an index which sets 
the efficiency, and /h 2 = Eh 2 / (Em + ^h 2 is the molecu- 
lar gas fraction. Table [1] gives a summary of the parameter 
values and their origin. 

Given a star-formation surface density, this equa- 
tion can be solved for the cold gas surface density when 
one knows the molecular fraction of the cold gas. We 
use a pressure-regulated recipe to determine the molec- 
ular fraction of the cold gas, based on the work by 
iBlitz fc Rosolowskvl (|2006l ). They found a power-law rela- 
tion between the mid-plane pressure acting on a galaxy disc 
and the ratio between molecular and atomic hydrogen, i.e., 




with Pq the external pressure in the interstellar medium 
where the molecular fraction is unity, a and fcs are the 
power-law index and the Boltzmann constant, respectively 
(see Table [T]). P m is the mi d-plane pressure acting on the 
galaxy disc, and is given by jElmegreenl ll98^) 



P m (r) = - GE gas (r) (E gas (r) + / CT (r)E„(r)) (3) 

where G is the gravitational constant, r is the radius from 
the galaxy centre, and / CT (r) is the ratio between <7 gas (r) 
and cr*(r), the gas and stellar vertical velocity di spersion, 
respec tively. We followed the method presented in lFu et al.l 
(2010) to determine f a as a function of radius r. Observa- 
tions have shown that cr gas remain s constant over th e disc 
(e.g. lShostak fc van der KruitJll984l ; lLerov et al.ll200Sl ). The 
stellar vertical velocity dispersion cr, decreases exponentially 
and has a scale len gth twice the stellar disc scale length r* 
(e.g. lBottemalll993l ). i.e. 

a*(r) = a°exp(-r/2r») (4) 
where the stellar velocity dispersion in the center of the disc 



is noted as ct*. The stellar mass surface density of an expo- 
nential disc is given by E* = E|J exp (— r/r*), which we can 
combine with the previous equation to find 

'-M = ^^ s rf^ (5) 

where /° is the ratio between the gas and stellar vertical 
velocity dispersion at the centre of the disc and E^ is the 
stellar mass surface density in the centre of the disc. We can 
express /° in terms of the mean value of f a (r) in the disc 
by 

- j27rrE.(r)/ g (r)dr 

U J 27rrE* (r)dr 1 ' 

Substituting for /o-(r) and afterwards integrating we find 

U = 4/°. (7) 
By combining Equation (0 and Equation {(Sj we find 

P m W = ^GE gas (r)(E gas (r) + ^y / ETwE2) (8) 

We now have all the necessary ingredients to calculate Ru 2 
and subsequently the cold gas molecular fraction [/h 2 = 
RkJ(1 + Rk 2 )]. 

Putting this together, we have an expression for the 
star formation law depending on the cold gas surface den- 
sity, stellar mass surface density and stellar mass and scale 
length of the disc. Given a star formation surface density 
and stellar surface-density, we can now indirectly determine 
the cold gas surface density through iteration. The molecu- 
lar gas fraction, and subsequently the molecular gas surface 
density, can then easily be calculated from Equation @, 
using the stellar mass surface density and previously deter- 
mined gas surface density as input. 

Intrinsically hidden in this method is the CO-to-H2 con- 
version factor (Xco)- This conversion factor was used to 
determine 1) the normalization constan t Asf for the em - 
pirical star formation law based on H 2 (|Bigiel et al.ll2008l ) 
and 2) the ratio between molecular and atomic hydrogen 
Ru 2 which was used to determine the normalization and 
power-law index fo r the press ure-based H2 fraction recipe 
l|Lerov et al.l |2008| ). However, IBlitz fc Rosolowskvl (|2006l ) 
proposed that changes in the CO-to-tb conversion fac- 
tor are unlikely to substantially alter results for the lat- 
ter point. Unless stated differently, we adopted a value of 
Xco = 2.0 x 10 20 cm- 2 (Kkms" 1 )" 1 . 
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Figure X. Comparison between indirect and direct measures of 
molecular gas surface densities, taken fro m the radial profiles of 
the THINGS galaxies iLerov et alj|2008h using our method. All 
datapoints are coded by galaxy and represent surface densities 
sampled along their host galaxy radius. 

2.2 Constructing a galaxy 

The method described in the previous subsection allows one 
to calculate the gas surface density, given a measured stellar- 
mass surface density and SFR surface density. However, for 
many high-redshift objects only global parameters are avail- 
able, combined with some morphological information (e.g., 
stellar scale length r*, half light radius fo.s). 

In order to calculate the cold and molecular gas con- 
tent of these galaxies, we distribute the stellar mass into an 
exponential disc following E*(r) = S2exp(— r/r*) where r* 
is the scale length of the stellar disc and E* is the stellar 
mass surface density in the centre of the disc (defined as 
E$ = M*/(2irrt) where M„ is the galaxy stellar mass). The 
star formation in the galaxy is distributed similarly, with a 
disc scale length 

r gas — Xgas Vtt , (9) 

where Xgas is a scale radius of the gas disc relative to the 
stellar disc (see Table [T]). The use of Xgas is motivated by the 
radia l extent of Hi gas in spiral galaxies (Br oeils fc Rhed 
1997). A galaxy's total cold gas budget can be found by 
integrating over the inferred gas surface densities. We find 
the best fit results for integrated masses (see Section 13. 2\ 
when distributing the stellar mass and star formation over 
an exponential disc out to five times the stellar disc scale 
length. 



3 CALIBRATION OF THE METHOD 

In this section we compare indirect measures of gas surface 
densities and integrated masses obtained using the method 
presented in the previous section to direct measures from 
the literature. Results were obtained using the parameter 
values as summarized in Table [1] Note that the purpose of 
this method was to measure the global cold gas and molec- 
ular content of galaxies, rather than local surface densities. 
Parameters were chosen to best achieve this purpose. We 



discuss the effects of changing the parameters in Section 

El 

3.1 Local galaxy surface densities 

Figures[T]and[2]show total cold gas and molecular gas surface 
densities obtained using our method. The indirect measures 
were computed from the radial profiles of the star- formation 
surface densities and stellar mass surface densities, in com- 
bination with the host g alaxy total mass a nd stellar scale 
length all obtained from ILerov et all (|2008l ). These predic- 
tions are plotted against literature values for t otal cold gas 
and m olecular gas surface densities taken from ILerov et all 
|200Sh 

We find good general agreement between results ob- 
tained using our method and literature values at molec- 
ular surface densities ranging from £h 2 ~ 1 to Eh 2 ~ 
300 Mq pc~ 2 (Figure [1]). There is a scatter around the one- 
to-one line of 0.1 to 0.3 dex at the lowest and more inter- 
mediate surface densities, respectively. We overestimate the 
lowest molecular surface densities for NGC 6946 and NGC 
5194, whereas we under-predict at high surface densities for 
these galaxies (also see NGC 3627). The star- formation ef- 
ficiency (SFE, the ratio of SFR over molecular gas) of these 
galaxies decreases significantly towards smaller radii (at the 
highest surface densities) compared to the predicted SFE 
fr om the mid-plane pressure on the disc (see Appendix F 
of ILerov et~alf 120081 ). The low SFE of these galaxies with 
respect to theoretical values will result in the observed sur- 
plus of molecular gas as predicted by our method. Opposite 
behavior is observed for these galaxies towards larger radii, 
where the molecular surface density is much lower. In this 
case, deviations between observed and theoretical SFE re- 
sult in a deficit of molecular hydrogen as computed by our 
method. 

We also find very good agreement between the total 
gas surface densities computed using our method and sur- 
face densities obtained from the literature (left-hand panel 
of Figure with smaller scatter especially at the lower 
gas surface densities. We predict a deficit of total cold gas 
compared to literature values at the highest surface densi- 
ties. Similar to the molecular gas surface densities, we as- 
cribe this deficit to differences between theoretically pre- 
dicted and observed SFE. Overall we find a mean ratio of 
Egasindircct/Egas observed = 0.94 with a standard deviation of 
0.021 for our method. 

We plot indirect total cold gas surface densities ob- 
tained with the KS law as a function of the literature values 
in the right-hand panel of Figure[2]for comparison. Although 
the scatter around the one-to-one line of the two methods is 
similar, there is a clear offset from this line when using the 
KS law, resulting in a general deficit of gas over the entire 
range of total gas surface densities. We find a mean ratio 
E gas indirect /E gas observed = 0.79 with a standard deviation of 
0.019 when using the KS law for the indirect gas measure. 
This was expected fo r the lowest surfac e densities where the 
KS law breaks down (|Bigiel et al.l l2008l. which also appears 
in Figure [2] However, the offset is surprising at the higher 
gas surface densities at which the KS law is defined to re- 
late SFR surface densities with the total gas surface density. 
Putting the observed gas surface densities in Esfr — E gas 
space and comparing with the KS law shows that the data lie 
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Figu re 2. Indirec t versus direct measures of total cold gas surface densities, taken from the radial profiles of the THINGS galaxies 
llLerov et alj|2008h . Left and right panel using our and the KS indirect gas measure method, respectively. We find better agreement with 
literature values than the indirect KS method, especially at low surface densities. 
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Figure 3. Indirect vs. direct measures of galaxy m olecular gas 
mass es using ou r indirect gas method. Both low- |Lerov et alj 
l2008t) and high- jTacconi etal]|20ld; [baddi et alj|2010h redshift 
galaxies are displayed. Apart from four local dwarf galaxies (see 
Sec. 13.2b . we find good agreement between literature values and 
our method at all redshifts and molecular gas masses. 



slightly below the KS law. Nevertheless, this can not account 
for the breakdown at low gas surface densities. We therefore 
find that our method yields the biggest improvements at the 
lowest gas surface densities (i.e., the outer parts of galaxy 
discs) . 



3.2 Galaxy integrated gas masses 

We present galaxy integrated Eh and total cold gas masses 
obtained with our method versus literature values in Figures 
[3] and 31 respectively. We applied our method to a small 
sample of bo t h low and high- redshift galax i es, com piled from 
iLerov et all (120081 . z = 0), iDaddi et all l|20ld . z ~ 1.5), 
and iTacconi et al.l <|201Ct z ~ 1.2 and z ~ 2.2). We applied 
the Xco factors as determined by the latter two authors 
for their galaxy samples [Xco ~ 2.5 x 10 20 and Xco = 
2.0 x 10 20 cm" 2 (Kkms -1 )" 1 , respectively]. 

We find good agreement between our results and ob- 
servations over the entire range of H2 masses and redshifts. 
There is a slight excess of H2 computed with our method as 
compared to the direct measurements (~ 0.2 dex; Figure |3J). 
We find a spread of ~ 0.25 dex around the mean trend in 
computed versus literature H2 masses. We neglect the four 
local galaxies with a surplus of integrated H2 masses of 1- 
1.5 dex. All four are classified as dwarf galaxi es and have 
anom alously high SFE in their central parts (|Lerov et alj 
2008). These authors argue for an unaccounted-for reservoir 
of H2 in these galaxies, possibly due to variations in the 
Xco fac tor for dwarfs with respe ct to bigger spiral galaxies 
(see also lSpaans fc Norman|[l997f ). We find one clear outlier 
at higher redshift, with an upper limit on its observed H2 
mass. It is most likely that this galaxy also has a high SFE 
we cannot account for. 

Current instruments do not allow for atomic hydrogen 
observations at high redshifts. We are therefore restricted to 
local galaxies when comparing literature values to total in- 
tegrated galaxy gas masses computed with the method pre- 
sented in this study. We find good agreement between our 
results and the observational data, although with a slight 
mean deficit (left-hand panel of Figure |4}. We found a max- 
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imum deviation in literature values and our results of ~ 0.5 
dex. Predicted cold gas masses for most galaxies are within 
~ 0.3 dex of their literature values. It is worth noting that 
the four dwarfs for which we over-predict integrated H2 
masses are also included in this figure and have predicted 
total gas masses in good agreement with the observations. 
Although this seems promising, it means that not only do we 
find a surplus of H2 compared to current literature studies, 
we also find a deficit of atomic hydrogen. 

In the right hand panel of Figure [4] we present indirect 
total cold gas measures for local galaxies obtained by apply- 
ing the KS law. Similar to the surface density profiles, the 
KS law predicts a deficit of total integrated gas mass. The 
outer parts of the galaxies (where the gas surface densities 
are lowest, the KS law breaks down and more gas is neces- 
sary to sustain SF) drive this deficit. Here we clearly see the 
strength of our method over previous indirect gas methods. 

These results demonstrate out that our method can be 
used as a useful tool to predict the cold and molecular gas 
surface densities, as well as integrated gas masses for low- 
and high-redshift galaxies. 

3.3 Changing parameters 

The method presented in this paper has several free param- 
eters, which all individually affect the outcome of indirect 
gas measures. In this subsection we explore the spread in 
results due to variations in free parameters. 

We find that our method is robust against changes in 
the mean ratio between gas and stellar vertical velocity dis- 
persion fa. We explore values ranging from 0.1 to 2.0 and 
find differences in total cold gas surface density measures 
up to ~ 0.1 dex with respect to results obtained using the 
quoted value for f a in Table (JXJ) . Differences are most promi- 
nent at the lowest gas surface densities (S gaa ^ 10M© pc -2 ). 
The spread in integrated total cold gas masses due to 
changes in / CT also reaches only ~0.1 dex over the entire mass 
range probed by our control sample of galaxies. Changing f a 
does not significantly affect molecular gas surface densities 
and integrated masses. 

iBlitz fc Rosolowskvl (I2006T I derived a power- law index 
a and a normalization constant Po slightly different than 
those we use (a = 0.8 versus 0.92 and P = 5.93 x 10~ 13 
versus 2.35 x 10 -13 , respectively). Applying their values to 
our method results in small changes of only ~ 0.1 dex at 
lower gas surface densities with respect to indirect gas mea- 
sures obtained using the parameters as summarized in Ta- 
ble [T] Integrated gas masses differ within ~ 0.15 dex over 
the entire mass range probed. 

The most significant changes in inferred indirect gas 
measures are due to variations in the n ormalization of the 
star formation law Asf - In recent works iBiriel et all C2008) 
and iBigiel et al.l (|2011h found values for the star forma- 
tion normalization of A$f = 8.42 x 1O _3 M yr -1 kpc -2 and 
Asf = 4.6 x 1O _3 M0 yr -1 kpc -2 , respectively (after con- 
version to a Chabrier IMF). We explored a range of pa- 
rameters around these values and applied a x 2 test on the 
derived integrated molecular and atomic gas masses (ex- 
cluding dwarf galaxies). This resulted in a best fit value 
of A S f = 8.4 x 1O" 3 M yr" 1 kpc~ 2 for the integrated to- 
tal cold gas masses, and Asf = 13.5 x 1O _3 M0 yr -1 kpc -2 
for the integrated molecular gas masses. We found that 



Asf = 10.9 x 10 _3 Mq yr -1 kpc -2 leads to best overall agree- 
ment between direct and indirect measures for both the total 
cold and molecular gas integrated masses (see Table [TJ. 

Adopting other values (running from Asf = 1.0 x 10 -3 
to Asf = 9.0 x 10 _4 M Q yr -1 kpc -2 ) for the normaliza- 
tion of the star formation law leads to differences of a 
few tenths of a dex between literature values and our to- 
tal cold gas and molecular gas results. For ex ample, Asf = 
4.6 x 10 _3 A/ Q yr _1 kpc" 2 jBigiel et al.ll201ll) results in to- 
tal gas surface densities differing by ~ 0.2 dex from results 
obtained using the parameter values presented in Table Q] 
This is a systematic shift upwards in indirect surface density 
over the entire range of surface densities probed. Changes 
in integrated indirect mass measures are of the order ~ 0.15 
dex with respect to results obtained using the tabulated pa- 
rameter values. 

Offsets between molecular gas direct and indirect mea- 
sures are more prominent, up to 0.5 dex for the inte- 
grated molecular gas masses when we adopt Asf = 4.6 x 
10 _3 Mq yr -1 kpc -2 . Adopting Asf with values larger than 
presented in Table [1] results in a similar systematic shift 
downwards for predicted surface densities and integrated 
masses. 

We applied an adapted formalism of the star forma- 
tion law presented in lBigiel et al.l ()2008f ) to allow for higher 
SFE in dense re gions (motivated by iDaddi et al.l |2010| ; 
iGenzel et al.|[2O10h . Adopting instead the standard SFE has 
no significant influence on the integrated cold gas and molec- 
ular gas predictions of local galaxies. Less concordance is 
reached between literature values and our predictions for 
the high redshift objects. We find a shift of approximately 
~ 0.2 dex upwards away from the one-to-one line with re- 
spect to results obtained with the inclusion of an increased 
SFE at highest gas surface densities. 

One of the key assumptions in this method is an expo- 
nential distribution of matter in the galaxy discs. SF occurs 
in molecular clouds (local clumps in the disc not following 
an exponential distribution) , which could lead to a local un- 
derestimation of the cold gas surface density. Nevertheless, 
an exponential disc seems a valid app roximation for star - 
forming galaxies on the main sequence l|Wuvts et al.ll201ll ). 
and the local clumps should average out when integrating 
over the disc. As we will discuss in Section [4.11 part of our 
COSMOS galaxy sample consists of quenched galaxies which 
fall off the main sequence. We are aware that these galax- 
ies are not described by disks but that their light profiles 
are b etter described by a Vaucouleurs profile l|Wuvts et al.l 

We assume that the gas in galactic discs is radially more 
extended than the stars. Varying the scale length of gas rel- 
ative to stars, Xgas ^ 1, results in only subtle changes in the 
integrated total cold and molecular gas masses. However, de- 
creasing the gas scale length to Xgas 

< 1 lowers the inferred 
integrated cold (molecular) gas masses. This difference can 
increase to ~ 0.25 (0.1) dex when Xgas = 0.5 with respect to 
results obtained using the parameters presented in Table [T] 
Both the SFR and stellar mass are a function of the 
IMF. Changes in inferred total cold gas surface densities are 
less than 0.1 dex at the lo west gas surfa ce densities when 
adopting a Salpeter IMF (|Salpeterl Il955h rather than the 
Chabrier IMF used above, and these changes become neg- 
ligible at the highest surface densities. The integrated cold 
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Figure 4. Indirect vs. direct measures of galaxy cold gas masses using our (left panel) and the KS (right panel) indirect gas method. 
No neutral gas measurements are available for the high-redshift galaxies, therefore only low-redshift galaxies are plotted. Our method 
finds good agreement with literature values, whereas the inverse KS method find a deficit of total cold gas mass. 



gas mass of the galaxies, as well as molecular gas surface 
densities and integrated masses do not change significantly 
when adopting a Salpeter IMF. 



4 GAS EVOLUTION IN THE COSMOS 
SAMPLE 

4.1 Sample 

In order to determine the cold gas content and atomic and 
molecular composition of high-redshift galaxies, we used 
a galaxy sa mple taken from the Cosmic Evolution Survey 
{COSMOS ; IScoville et al1l2007l ). The survey was designed 
to probe the evolution of galaxies, star formation, active 
galactic nuclei, and dark matter over the redshift range 
z > 0.5 — 6. These data comprise imaging and spectroscopy 
from X-ray to radio wavelengths including HST imaging and 
cover a two-square degree area on the sky. 

We considered a catalogue of ~ 57000 galaxies with 
HST / -band magnitude Iab < 24 mag in the redshift range 
0.5 < Zphot < 2.0 from the COSMOS survey. P hotometric 
redshifts were obtained from lllbert et all i|2009h . SFRs are 
based on rest-frame UV fluxe s, corrected using reddening 
values E(B-V) obtained from lllbert et all (|2009h for each 
indivi dual galaxy . UV fl uxes were converted into SFR fol- 
lowing [Kemficutt] fl998a) for a Chabrier-IMF. We computed 
stellar masses using a multi- wavelength SED \ 2 fitti ng to the 
COSM OS photometry (Uto _R"-band), applying the iBruzuall 
templates fixed at the photometric redshifts of the 
sam ple objects. We used g alax y half light radii (r so) taken 
from IScarlata et al.l (|2007T ) and lTasca et al.l (|2009l ). 

The selection of galaxies in the /-band does not nec- 
essarily guarantee that all selected galaxies lie on the M*- 
SFR plane. Our selection rather consists of both active star 
forming galaxies and galaxies already quenched. Inspection 



of the data points out that quenched galaxies indeed are 
especially dominant towards the massive end of our galaxy 
sample (logM* > 11 [M ]). 

We warn the reader that the sampling of galaxies at 
2 > 1.5 (and especially z > 1.75) is sparse, particularly at 
the low-mass end and one should use caution interpreting 
results in this redshift range. 



4.2 Cold gas content 

We present the derived total cold gas and molecular gas 
masses in Figures [5] and [6] Grey scales represent the log 
of the number of galaxies in each bin in the M gas {Mh2)~ 
M» plane, with 50, 16 and 84 percentiles over-plotted on 
the figure. Results are presented as a function of stellar 
mass in different redshift bins. Cold gas masses range from 
log Afgas ~ 10-10.5 Mq at z = 0.5-0.75 with a peak around 
logM* ~ 9.5 M e to log Af gas ~ 11-12 M Q at the high- 
est redshifts in our sample peaking at a stellar mass of 
logAf, ~ 10.3 M Q . 

We find a large spread in molecular gas masses in the 
range of logA/f/2 ~ 8.0-11.0 Mq at the lowest redshifts in 
our sample and logAfn2 ~ 10.0-11.5 Mq for the highest- 
redshift objects in our sample. The molecular gas masses 
peak at the same stellar masses as the integrated total cold 
gas masses. 

Not only do the cold and molecular gas masses of galax- 
ies at fixed stellar mass decrease with time, the largest gas 
reservoirs are always found in less- (stellar) massive objects. 
This trend suggests that (stellar) massive galaxies consumed 
their cold gas first, and only at later stages do less- massive 
galaxies fully consume their gas reservoir. We have explored 
the SFRs of galaxies of our sample galaxies as a function 
of redshift and found that for a given stellar mass, higher- 
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Figure 5. Total-cold mass gas as a function of stellar mass for different redshift bins. The grey shaded area shows the log of the number 
of galaxies in each gas/stellar mass bin, with the 50, 16 and 84 percentile curves shown with the red solid and dashed lines. The central 
bottom panel shows the 50 percentile curve for the data in each redshift bin. Bumps at the high-mass end in the highest redshift bins 
are due to ULIRGS (see Sec. IP) 
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Figure 6. H2 mass as a function of stellar mass for different redshift bins. The grey shaded area shows the log of the number of galaxies 
in each bin, with the 50, 16 and 84 percentile curves shown with the red solid and dashed lines. The central b ottom panel shows th e 50 
percentile curve fo r the data in each redshift bin. Red squares and pentagons are literature values taken from lTacconi et al.l j2010l) and 
iDaddi et al] I^OlOh . respectively. 
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Figure 7. Total cold gas fraction (atomic plus molecular) as a function of stellar mass for different redshift bins. The grey shaded area 
shows the log of the number of galaxies in each bin, with the 50, 16 and 84 percentile curves shown with the red solid and dashed lines. 
The central bottom panel shows the 50 percentile curve for the data in each redshift bin. The dashed and dotted yellow lines represents 
the fit and one sigma scatter of Equation 1101 respectively. 



redshift objects have higher SFR (i.e., consume their gas 
faster) than their low-redshift counterparts. 

We have superim pose d observational result s obtained 
by iDaddi et all l|201Ch and iTacconi et all l|201Cj| on Figure 
[U We obtain good agreement between molecular gas masses 
computed with our method and observations extracted from 
the literature. Only in the redshift regime z = 1.0-1.25 
do we compu t e a d eficit of Hb compared with results of 
ITacconi et ail |2010l ). 

As pointed out in section 14.11 quenched galaxies make 
up for the dominant part of the massive obje cts in our sam- 
ple an d do not track the M*-SFR plane. The ITacconi et all 
(|2010f ) galaxies on the other hand are actively forming stars 
and lie on the M*-SFR plane. We find that these litera- 
ture galaxies have specific SFR (SSFR) much higher than 
galaxies in the same redshift bin and with similar stellar 
masses in our COSMOS sample (> 1.0 dex). This difference 
in galaxy population (and consequently SSFR) can explain 
for the surplus of Eb in literature galaxies compared to our 
predictions. 

We find a few massive objects (logM* ~ 11.5— 
12.0 Mq) with large total cold and molecular gas reservoirs 
(logMg as ~ 12 Mq), in disagreement with the declining 
trend in cold gas mass in less (stellar) massive objects. These 
are a handful of objects with SFRs of up to ~ 2OOM0yr _1 
and drive the increment of the 50 percentiles in the highest 
redshift bins. The objects are ultra luminous infrared galax- 
ies (ULIRGS) with bright 24 /im intensities (~ 0.5 mjy), 
and are likely to have an AGN component. There are more 
ULIRGS in our sample, but it is the combination of high 
stellar mass and SFR which makes these galaxies show up so 



prominently. If these galaxies indeed have an obscured AGN 
component, both the stellar mass and SFR could be overes- 
timated which would be refle cted in the cold an d molecular 
gas content of these objects (|Caputi et al.1 12006). 



4.3 Gas fraction 

We present the gas fraction (fraction of all the baryonic mass 
which is in cold gas) of our galaxy sample as a function of 
stellar mass for different redshift bins in Figure [7] We find 
a clear trend in gas fraction and stellar mass for all redshift 
bins, with the highest gas fractions at low stellar masses. The 
gas fraction drops towards higher stellar masses and remains 
fairly constant around 0.0-0.1. More prominent, however, 
is the evolution of this trend with redshift (see the central 
panel in the third row of Figure[jJ). The stellar mass at which 
the gas fraction drops increases with increasing redshift, or, 
in other words, the gas fraction of galaxies at a given stellar 
mass increases with increasing redshift. 

Figure [8] shows the fraction of molecular hydrogen H2 
with respect to the total stellar and molecular gas mass 
(Mh 2 /(Mu 2 +M*)) as a function of stellar mass. This allows 
for a direct comparison with high-redshift galaxies from the 
literature. We find a trend between Mh 2 / (Mh 2 + M*) and 
stellar mass similar to the results of Figure [7J It also shows 
that A1h 2 /(Mh 2 + M») at a fixed stellar mass increases with 
redshift. There is a large spread in Mh 2 /(Mh 2 +M*) at fixed 
stellar mass, especially in the lower redshift bins. As in Fig- 
ure [6l we find good agreement between observations in the 
redshift range 1.25 < z < 1.75 and our results. We find a 
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Figure 8. A/h 2 /(Mh, + M,) as a function of stellar mass for different redshift bins. The grey shaded area shows the log of the number 
of galaxies in each bin, with the 50, 16 and 84 percentile curves shown with the red solid and dashed lines. The central bottom panel 
shows the 50 percentile c urve f or the data in each redshift bin. Red squares and pentagons are literature values taken from lTacconi et alj 
j201Ch and lDaddi et al.1 1 I2OIOI ). respectively. The dashed and dotted yellow lines represents the fit and scatter (one sigma) of Equation 
1111 respectively. 
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deficit in Mh 2 /(Mh 2 + M„) compared to observed objects at 
lower redshifts (1 < z < 1.25). As explained in the previous 
subsection this can be caused by differences in galaxy popu- 
lation. The peak in cold gas fraction and A1h 2 /(Mh 2 + M*) 
in the two highest redshift bins at high stellar masses rep- 
resent the few galaxies with high SSFR which are classified 
as ULIRGS. 

The cold gas fraction in our sampled redshift range can 
be characterized by a function of stellar mass and redshift, 
i.e., 



Mgas + M, exp( lo s M * ~ A 1/ B +1 ' 



(10) 



where A = 9.03(1 + q%)° 11 and B = 0.42(1 + z)" ' 56 . 
This fitting formula was designed to match the 50 per- 
centiles in 50 redshifts bins between 0.5 < z < 1.7 (where 
we sampled a larg e nu mber of galaxies) and lo cal data from 
iLerov et al.l (|2008h and lSaintonge et ail (120111 ). respectively. 
We included local data with the aim of producing a rep- 
resentation of galaxy gas fraction valid for both local as 
distant galaxies. We find an approximately constant scat- 
ter of o = 0.18. This slightly decreases towards the highest 
redshifts, but is likely an underestimation due to the low 
number of massive objects with large variations in their gas 
fractions. Although the function was not designed to match 
the highest redshift objects, it does a good job in reproduc- 
ing their gas fraction. The left-hand panel of Fig.[9]shows the 
distribution of gas fraction with mass of local galaxies. We 
superimpose our fitting function and find that it describes 
the observed data fairly well. 

Using the same approach we designed a fitting function 
to Mh 2 /(Mh 2 + M t ), which is given by a similar equation 



M H , 



Mh 2 + M, expC lo s m*-a)/b +1 ' 



(11) 



where A = 6.15(1 + o^) ' 144 and B = 1.47(1 + z)~ 2 - 23 . 
This formula is a good quantitative representation of our 
full sample for local galaxies up to z — 2, and has a scatter 
of a = 0.21 for our high-redshift galaxy sample (much less 
for the local objects). 

Note that at stellar masses above logA/„ > 11.0 [M©] 
the presented fitting formulae do not apply for active star 
forming galaxies, but rather for an 'overall' galaxy pop- 
ulation which consists of a mixture of star forming and 
quenched galaxies. We advise the reader to use the method 
presented in section [5] in the regime of massive, actively star 
forming galaxies. 



redshifts [z = 1-1.5) and remains at this level at higher 
redshifts. The most-massive objects at the highest redshifts 
(z — 1.5-2.0) have molecular gas fractions of fn 2 = 0.2— 
0.3 and, similar to the galaxy gas fraction, show a small 
peak at log M* /Mq ~ 11.5-12.0 (caused by the ULIRGS). 
In general, massive galaxies reach low /h 2 values before less 
massive galaxies do. This trend is connected to the evolution 
of the gas fraction discussed in the previous subsection. 

The described trend is reversed compared to the trend 
observed for local galax ies (i.e., decreases with de- 

creas ing stellar mass, e.g.. lLerov et al.ll200Sl : ISaintonge et al.l 
l2011f ). This disagreement is driven by the quenched galax- 
ies (with relatively low SFR and consequently relative low 
molecular content) which are dominant towards the mas- 
sive end of our galaxy sample. We can decouple the mixture 
of passive and actively star-forming galaxies by considering 
galaxies at fixed SSFR. We should then expect to see the 
trend in fn 2 with stellar mass found in local galaxies, pre- 
cisely as seen in Figure [TT1 

The molecular fraction of the cold gas with respect to 
the total cold gas in the galaxies for different redshift bins 
is presented in Figure 1121 The figure shows a peak in H2 
fractions at higher redshift. At the lowest redshifts we see a 
gentle increase of H2 fraction with increasing total cold gas 
mass, whereas at higher redshift the H2 fraction peaks at 
a total cold gas mass of logM gas ~ 10.7 and decreases at 
higher masses. 

This behavior demonstrates that the apportioning of 
the cold gas into atomic and molecular hydrogen is not 
driven by the cold gas budget availab le. This is not surpris- 
ing as the lBlitz fc Rosolowskvl | |2006t ) formalism is based on 
the cold gas surface density. 



4.5 Transition between gas and stellar dominated 
phase 

The transition from a gas-dominated phase to a stellar- 
dominated phase can be used as a probe of the evolution- 
ary stage of a galaxy. We plot the stellar masses of galaxies 
with gas fractions around ~ 0.5 (which marks the transition 
from gas to stellar dominated) as a function of redshift in 
Figure [131 We find that the typical stellar mass at this tran- 
sition point decreases from log M» (f glLS /MQ = 0.5) = 11.1 
at * = 2.0 to log M„ (/gas/M© = 0.5) = 10.1 at z = 0.5. This 
decrease of an order of magnitude in stellar mass indicates 
that massive galaxies reach the transition from a gas- to a 
stellar-dominated phase well before less-massive galaxies. 



4.4 H2 fraction 

Figure [TU] shows the H2 fraction of the cold gas (fraction of 
cold gas which is in a molecular state, H2/(H2 +HI)) in our 
galaxy sample as a function of galaxy stellar mass, divided in 
different redshift bins. Galaxies have H2 fractions spanning 
a range from nearly zero (hardly any H2) to almost one (all 
the gas is in molecules) at all redshifts. There is nevertheless 
a (median) trend for the overall population of the galaxies. 

We find the molecular fraction of the gas in most galax- 
ies up to z = 1 to be around fn 2 = 0.2, with a minor de- 
crease of the H2 fraction with increasing stellar mass. The 
H2 fraction increases up to values of fn 2 = 0.3-0.4 for most 
low-mass objects (below log M,/Mq ~ 10) at intermediate 



5 DISCUSSION 

We have shown that one can fairly accurately predict the 
total cold gas content of star forming galaxies as a function 
of their SFR, stellar mass and size. The cold gas can fur- 
ther be apportioned into atomic and molecular hydrogen. 
The key ingredient is a molecular gas-based star formation 
law combined with a prescription for the calculation of cold 
gas molecular fractions. In order to model a galaxy, the 
method assumes an exponential disc, integrated outwards 
to five times the stellar disc scale length. We obtain the 
best results when increasing the SFE in high cold gas sur- 
face density regions, especially important for high-redshift 
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Figure 10. The fraction of H2 in cold gas as a function of stellar mass for different redshift bins. The grey shaded area shows the log 
of the number of galaxies in each bin, with the 50, 16 and 84 percentile curves shown with the red solid and dashed lines. The central 
bottom panel shows the 50 percentile curve for the data in each redshift bin. 
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Figure 11. The fraction of H2 in cold gas as a function of stellar mass for galaxies with a fixed SSFR. All galaxies have redshifts 
0.75 < z < 1.0. The grey shaded area shows the log of the number of galaxies in each bin, with the 50, 16 and 84 percentile curves 
shown with the red solid and dashed lines. Left panel: galaxies with SSFR —9 < log SSFR < —8 [yr -1 ]. Right panel: galaxies with SSFR 
— 10 < log SSFR < —9 [yr -1 ]. Note that fu 2 decreases with decreasing stellar mass. 



galaxies (see Sec. 13. 3[) . This is in agreement with previous 
studies suggest ing an increased SFE for dense high-redshift 
objects ijDaddi et al.ll2010l : Icienzel et al.ll20TOt ). It is outside 
the scope of this work to address the physical origin of this 
change in SFE. 

Our method provides significantly better predictions of 
galaxy gas estimates than those derived from the KS law, 



as can be seen from the comparison with direct galaxy gas 
measurements extracted from the literature. This is due to 
the break-down of the KS law at low surface densities, re- 
sulting in an underestimation of the cold gas necessary to 
sustain low SFR-rates. This break-down leads to significant 
changes when integrating over the outer parts of exponential 
discs. 
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Figure 12. The fraction of H2 in cold gas as a function of total gas mass for different redshift bins. The grey shaded area shows the log 
of the number of galaxies in each bin, with the 50, 16 and 84 percentile curves shown with the red solid and dashed lines. The central 
bottom panel shows the 50 percentile curve for the data in each redshift bin. 
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Figure 13. Stellar masses of galaxies with / gas between 0.47 
and 0.53 as a function of redshift. Grey scales show the log of 
the number of galaxies in each bin. Massive galaxies reach the 
transition point from a gas to a stellar dominated phase before 
less massive objects do. 



Our method is applicable to low- and high-surface- 
density regions, as well as low- and high-redshift objects. 
Furthermore, the empirical success of the method gives sup- 
port to our assumption of exponential profiles for the stellar 
and gaseous discs. Conversely, the method developed seems 
to be a good approach to compute SFRs from gas and stellar 
masses. 

Despite the success of our method, it fails to properly 



predict the molecular and tot al cold gas c onten t of dwarf 
galaxies. As stated previously, iLerov et al.l (2008) argue for 
an unaccounted-for reservoir of H2 in these galaxies, possibly 
due to variations in the Xco factor for dwarfs with respect 
to bigger spiral galaxies. It would therefore be worthwhile 
to address the Xco conversion factor in dwarfs, in order 
to better understand the gas physics driving the partition 
of cold gas into atomic and molecular hydrogen. This will 
be crucial for a proper understanding of star formation at 
different stellar masses. 

It is worthwhile to remember that we assume that the 
Xco factor and IMF are universal. An evolving or variable 
IMF and/or Xco could potentially lead to significant biases 
in our estimates. A much lower Xco (suggested for starburst 
galaxies in e.g., iGenzel et al.|[201oT ) would result in much 
smaller H2 reservoirs than predicted by our method. The 
IMF enters in converting the observed luminosities to stellar 
masses, as well as converting the observational tracers (such 
as H-a, UV or IR luminosity) to star formation rate. A top- 
heavy IMF would result in higher UV, IR or H-a luminosities 
for a given SFR, or correspondingly lower SFR estimates 
for a fixed observed luminosity. This would result in lower 
derived gas fractions. 

Although the gas masses and molecular fractions pre- 
sented here are obtained indirectly, our method allows for 
the first time the quantification of the cold gas and molec- 
ular gas for a large sample of galaxies at high redshift. This 
allows us to study the evolution of gas content in galax- 
ies and track the relation between stellar and gas mass. 
High-redshift star-forming objects appear to be more gas- 
rich than lower redshift objects with similar stellar masses 
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(Fig. We found that objects at z ~ 2 have a cold gas 
reservoir approximately thirty times larger than objects at 
z ~ 0.5 (logAfgas/Af© ~ 11.5 versus logM g as/M ~ 10.0 at 
a stellar mass of log M*/Mq = 8.5) and that the largest gas 
reservoirs are found in progressively less-massive galaxies as 
a function of time, molecular gas reservoirs follow this be- 
havior closely. These results indicate that massive galaxies 
consume or expel their gas earlier and at a higher rate than 
less-massive galaxies. 

The galaxy gas fractions share the same evolution as 
total cold gas contents and decrease rapidly with decreasing 
redshift (Fig. [7|. This trend is in good agreement with the 
few direct high redshift molecular gas measures available to 
date (Fig. [8]). More striking is that massive galaxies reach 
low gas fractions (less than ~ 0.1) much sooner than less- 
massive objects. 

We find a weak trend between stellar mass and the frac- 
tion of gas in molecular form (see Figure [10)) . such that mas- 
sive galaxies have somewhat lower molecular fractions. In 
addition, there is a weak redshift trend in the sense that, at 
a fixed stellar mass, galaxies have higher molecular fractions 
at earlier times. 

The apparent weak trend between stellar mass and the 
fraction of gas in molecular form is driven by the mixture of 
galaxies on the main-sequence and already quenched galax- 
ies. We find that the molecular fraction of cold gas increases 
with stellar mass when we decouple this mixture (see Figure 
Hip. This is in agreement with observatio ns of local galaxies 
ijlTerov et al.ll2008l ; [Samtonge et al.ll201ll ). 

We find no clear relationship between the molecular gas 
fraction and the total mass of cold gas available (see Figure 
I12p . A large gas reservoir does not necessarily imply a high 
molecular fraction. However, we do find that the molecular 
fraction of a galaxy with fixed cold gas mass decreases with 
decreasing redshift. Under the assumption that cold gas is 
only eligible for star formation in a molecular state, these 
results indicate that high-redshift objects with a fixed cold 
gas content can transform more cold gas into stars than low 
redshift objects. 

As discussed in Section [3.31 and [ 4.11 part of our galaxy 
sample falls off the main-sequence and is best described by 
a Vaucouleurs light profile rather than an exponential disc. 
These objects push to the limits of our method (which as- 
sumes an exponential disc) and one should interpret these 
with care. 

In a hierarchical galaxy formation scenario, 
galaxy gas fractions represent the competition be- 
twee n gas inflow, outflow and consu mption through 
SF dDave. Finlator fc Qppenheimerl lioilT ). All the previ- 
ously discussed processes hint at a similar general behavior 
of this competition with respect to host galaxy redshift 
and stellar mass. Massive galaxies consume or expel their 
gas before less-massive galaxies do. At a fixed stellar mass 
galaxies consume or expel more gas than they attract. This 
is all in agreement with iDave. Finlator fc Qppenheimerl 
|201ll ). who argue that the cosmic inflow rate of gas 
diminishes faster than the consumption rate. Furthermore, 
we find that the molecular fraction of the cold gas decreases 
with both stellar mass and time. 

We present the transition from a gas to a stellar domi- 
nated galaxy in Figure [T3l Galaxies that pass this transition 
point have expelled or consumed most of their gas and are 



likely (unless they accrete a fresh supply of gas) to proceed 
towards a more quiescent phase of SF. We find that more 
massive galaxies reach a quiescent gas-poor state before less 
massive galaxies. 

These processes are a reflection of downsizing in star for- 
mation, which can be defined as the declini ng mass of star- 
form ing galaxies with decreasing redshift |Fontanot et al.l 
2009). Our results suggest that the slowly shrinking gas 
reservoirs, combined with lower molecular gas fractions, may 
drive the downsizing observed in star formation. It will be 
interesting to test this picture with future 'direct' observa- 
tions of gas in galaxies at high redshift. 



6 SUMMARY AND CONCLUSION 

We developed a method to indirectly measure the total cold 
gas and molecular gas content of galaxies. We applied this 
method to a sample of galaxies from the COSMOS survey 
covering a redshift range of 0.5 ^ z ^ 2.0. Since we obtained 
galaxy gas masses from their SFR, our results are effectively 
another way of representing SF and might seem redundant 
with other previous studies in the literature. It is therefore 
worthwhile to briefly summarize what we actually learned 
from this approach that we could not have concluded from 
the SFR alone. Our main results are: 

• The gas content of a galaxy can be described using a 
combination of a molecular SF-law and a prescription to 
calculate the molecular fraction of cold gas, all under the 
assumption of a radial exponential profile for the gas and 
stars. Best results are obtained using a variable SFE. Con- 
versely, this method can act as a scaling law relating the 
SFR surface density to the cold gas surface density of a 
galaxy. The method seems to be applicable to a large range 
of galaxy gas masses at both low and high redshift. 

• The method presented allows us not only to predict 
the cold gas content but also to predict the molecular gas 
content of galaxies, and it is based on observed galaxy prop- 
erties. These estimates can help us to interpret direct obser- 
vations of molecular gas in galaxies at high redshift relative 
to the overall galaxy population. 

• We find a strong trend between the gas content of a 
galaxy and its stellar mass with time. For a given stellar 
mass, the gas fraction of the galaxy decreases with decreas- 
ing redshift. On average, massive galaxies consume or expel 
their gas reservoir much earlier in the history of the Universe 
than less massive galaxies. This trend is not only observed in 
the cold gas content and gas fraction of the galaxies; it is also 
reflected in the molecular fraction of the cold gas, which gets 
smaller with time for a given stellar mass. This, in combi- 
nation with low gas fractions, shows that the physics which 
suppress the formation of stars in massive galaxies with time 
is at least two-fold: galaxies run out of gas and molecules, 
but not necessarily at the same rate. 

• The stellar mass at which galaxies become stellar- 
dominated decreases with time. This indicates that massive 
galaxies reach a gas-poor state before less massive objects. 

• These results point towards a common trend: the more 
massive a galaxy, the sooner it consumes or expels its gas 
content. This trend is another manifestation of downsizing. 



Gas evolution in galaxies at 0.5 ^ z ^ 2.0 15 



Although still indirect, these results point out that mea- 
suring the gas content of galaxies is of great use for devel- 
oping a broader picture and better quantification of galaxy 
evolution. The method and results presented in this work 
are a first step in this direction and can help for the de- 
sign of upcoming surveys. Furthermore, they provide useful 
constraints for cosmological galaxy formation models, par- 
ticularly for the development of new models which include 
detailed tracking of multi-phase gas. 
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